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Abstract — A recent way to model and analyze downlink cellular 
networks is by using random spatial models. Assuming user 
equipment (UE) distribution to be uniform, the analysis is 
performed at a typical UE located at the origin. While this 
method of sampling UEs provides statistics averaged over the 
UE locations, it is not possible to sample cell interior and cell 
edge UEs separately. This complicates the problem of analyzing 
deployment scenarios involving non-uniform distribution of UEs, 
especially when the locations of the UEs and the base stations 
(BSs) are dependent. To facilitate this separation, we propose a 
new tractable method of sampling UEs by conditionally thinning 
the BS point process and show that the resulting framework can 
be used as a tractable generative model to study cellular networks 
with non-uniform UE distribution. 

I. Introduction 

A promising new way to model cellular networks is by using 
random spatial models, where the BS locations are assumed 
to form a realization of some spatial point process, typically 
the Poisson Point Process (PPP). Modeling the UE locations 
as an independent PPP, the downlink analysis is performed 
at a typical UE assumed to be located at the origin. Owing 
to its tractability, this model leads to simple closed form 
expressions for key metrics such as coverage and average rate 
over the entire network Q]-|3). Although it provides a way 
to study average statistics over the complete network, it does 
not provide any handle on studying the performance of cell 
edge and cell interior UEs separately. Consequently, it is not 
possible to model deployment scenarios involving non-uniform 
UE distribution, especially when the UE and BS locations are 
dependent. Naturally, such flexibility is desirable, especially 
to study the current capacity-centric deployments where BSs 
are specifically deployed to be proximate to areas of high 
UE density. The most popular and perhaps the only available 
option to handle such scenarios is through detailed system 
level simulations, which are both time consuming and have to 
be focused on a limited range of system parameters [4|, [5|. 

As a first step towards a tractable model, we propose a slight 
modification in the way UEs are sampled by introducing the 
idea of conditional thinning. Starting with a higher density of 
BSs, we first assume that a typical UE is located at the origin. 
After selecting the serving BS, we condition on this active 
link and independently thin the rest of the BS point process 
so that the resulting density matches the desired density of the 
actual BSs. The thinning operation pushes the typical UE in 
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the cell interior relative to the new cell edge defined by the 
resulting point process. The bias induced in the location of a 
typical UE towards its serving BS can be tuned by varying the 
thinning probability. We make this notion precise by deriving 
the distribution of the ratio of the distances of the UE to its 
serving BS and the dominant interferer as a function of the 
thinning probability. We also show that this framework can be 
used as a tractable generative model to study non-uniform UE 
distributions where the UEs are more likely to lie closer to the 
BS. The exact analysis of such non-uniform UE distributions 
is in general hard due to the correlation present in the UE and 
BS locations. The impact of the proposed model on the cellular 
performance analysis and its key differences from the existing 
model based on uniform UE distribution are highlighted in 
terms of the coverage predictions. 

It is worth noting that although this work is developed in the 
context of cellular networks, it applies to much wider class of 
point process problems involving dependence in the location 
of the observation point and the point process. 

II. Proposed Method of Sampling UEs 

We consider a homogeneous PPP $ of density A, a thinned 
version of which will eventually model the BS locations. The 
downlink analysis is performed at a typical UE assumed to be 
at the origin (6j . The received power at a typical UE from the 
BS located at x e $ is 

P x = Ph x \\x\\~ a , (1) 

where P is the transmit power, h x ~ exp(l) models channel 
power distribution under Rayleigh fading and models 
standard distance based path loss with a > 2 being the path 
loss exponent. More general fading distributions can be studied 
using tools developed in |7), (8). For simplicity of exposition, 
we will ignore thermal noise in this discussion. 

The proposed method of sampling the interior UEs can be 
understood in two simple steps shown in Fig. [T] The first step 
is to identify the serving BS in <E> for a typical UE depending 
upon the cell association technique being considered. While 
the proposed method is general, for brevity we consider 
nearest-neighbor cell association model discussed in 1 1 1, where 
each UE connects to its nearest BS. It is also the same as 
maximum average power connectivity model, where each UE 
connects to the BS that provides maximum long-term average 
received power. In Fig. [T] (left), the typical UE connects to its 
nearest BS, which is the one that corresponds to the Voronoi 
cell in which it lies. Denote the location of this serving BS 
by x s £ $. The second step is to independently thin the point 




Fig. 1. (left) The Voronoi tessellation of the point process Dark triangle denotes a typical UE. (middle) The point process thinned by p 
remaining points form BS point process (right) The Voronoi tessellation of $5. Observe that the typical point is now in the cell interior. 
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process $ \ x s where each point is independently retained 
with probability p as shown in Fig. [T] (middle). Since the 
thinning is conditional on the serving BS, we call it conditional 
thinning. The thinned process <!>;, models the BS locations. 
Due to conditional thinning, the point process $b is no longer 
a PPP. After tessellating the space based on <!>(, and keeping 
the position of the UE fixed, we note that the UE is pushed 
towards the cell interior compared to the new cell edge as 
shown in Fig. [T] (right). As discussed in the next section, a 
typical UE can be pushed arbitrarily close to its serving BS 
by choosing an arbitrarily small value of p but remains edge 
biased for high values of p. 

ILL Impact on Cellular Networks 

After providing an overview of the new UE sampling model, 
we now formalize the effect of conditional thinning perceived 
by a typical UE. Denoting i?i and R 2 to be the distances of 
the closest and the next closest point of $b to the origin, we 
define our first metric to be R = R2/R1. It corresponds to 
the ratio of the distances of a typical UE to its serving BS 
and the dominant interferer, and provides some insights into 
the expected performance of a typical UE. For instance, if the 
value of R is close to 1, it means that the UE is near the cell 
edge, i.e., the dominant interferer is approximately at the same 
distance as the serving BS and hence the received signal-to- 
interference-ratio (SIR) is expected to be low. It is important 
to note that R > 1 by construction, since the serving BS is 
always the closest one in our cell association model. We now 
derive the distribution of R after deriving the joint distribution 
of i?i and R 2 as a function of thinning probability p in the 
following Lemma. Interested readers can refer to [9| for the 
marginal distribution of R n . 

Lemma 1. The cumulative distribution function ( CDF) of R 
is 



F R (r) = 1 - 



1 



r > 1. 



(2) 



1 +p(r 2 - 1) 

Proof: For notational simplicity, we first define the fol 
lowing disjoint sets: 

£i = {x <E R 2 : \\x\\ < n} 



£ 2 = {xe 

£i = {x<E 
£4 = {x E 



n < \\x\\ < r x + dri} 
n + dri < ||a;|| < r 2 } 
r 2 < \\x\\ < r 2 + dr 2 }, 



where £\ denotes a circle centered at origin and the rest denote 
annular regions defined by concentric circles centered at the 
origin. Further let N(£) be a random counting measure of a 
Borel set £, i.e., N(£) = # of points in £. Now to derive the 
CDF of R, we first derive the joint probability density function 
(PDF) of R\ and R 2 , which by definition can be expressed as 

' (Ri & £ 2l R 2 S £4) 



hm 

dri->0 
dr 2 ->0 



(3) 



dridr 2 

The numerator of the above expression can be expressed as: 

P(i?i e £ 2 ,R 2 e £4) = 

P [N{£\) = 0, N{£ 2 ) = 1, N(£ 3 ) = 0, N(S 4 ) = 1) = (4) 
P (N{£i) = 0) P (N(£ 2 ) = 1) P (N(£ 3 ) = 0) P (N(£ 4 ) = 1) , 

where the simplification follows from the fact that the sets £ j 
are disjoint. Now recall that the point process $ is indepen- 
dently thinned by probability p outside the circle of radius 
fj + dri. Therefore, the above expression can be written as: 

F(R X € £ 2 ,R 2 e £4) = e- x * r * 

\tt [(n + dn) 2 - r 2 } e -M(n+dn) 2 -,?] e - P A4. 2 Mr 1+ dn) 2 ] 
p\TT [(r 2 + dr 2 ) 2 - r 2 2 ] e -^[(r 2+ dr 3 f -rj] ? (5) 
which for vanishingly small dri and dr 2 can be simplified to: 

p(2n\) 2 rir 2 exp (— A7rr 2 (l — p)) exp [—pX-nr 2 ^ dr\dr 2 , 

(6) 

from which the joint PDF of i?i and R 2 can be expressed as: 

/ Rl ,fl 2 (ri,r 2 ) =p(27rX) 2 r 1 r 2 e- x ^ 1 -^e-^ r K (7) 

for r 2 > ri > 0. Using this joint density, the CCDF of R can 
now be expressed as: 



> r = 



R2 
Ri 



> r 



(8) 




Fig. 2. The CDF of R for various values of conditional thinning probability 
p. The theoretical results are overlaid with dotted plots of the simulation 
results showing perfect match. 
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(9) 



which completes the proof. ■ 
Recall that for p = 1 the proposed model reduces to the 
uniform distribution of UEs, leading to the following corollary. 

Corollary 1. The CDF of R for the typical observation point 
that is defined to be uniformly distributed in M 2 independent 
of the BS locations (p = 1 in the proposed model) is 

1 



F R (r) = 1 



r > 1. 



(10) 



Corollary 2. The mean value of R as a function of p is 
1 ^-taa-^-V? 



E[R] = 1 



which for p = 1 is E[R] = 2 



VT^p 



(11) 



Remark 1 (Dependence in UE and BS point processes). The 
level of dependence induced in the locations of the UEs and 
the BSs is inversely proportional to the thinning probability p, 
i.e., the probability of finding a typical UE in the cell interior 
close to its serving BS is higher for smaller values of p. This 
is evident from the mean of R given by Corollary^and from 
the plot of the CDFs of R for various values of p in Fig. [2] 

We now derive the coverage probability of a typical UE 
sampled under the new proposed method. The coverage prob- 
ability P c denotes the average fraction of UEs in coverage and 
can be formally defined as the CCDF of SIR as follows: 



1 (SIR > /?) 



h Xa \\x s \\~ 



E 



ve* b \x. MVll 



>0\, d2) 



where the serving BS is assumed to be located at x s . Using 
tools developed in [1|, a simple expression can be derived for 
the coverage probability of a typical UE under the proposed 
sampling method. The main result is given in the following 
Lemma along with a brief proof sketch. 



Lemma 2. The coverage probability of a typical UE when it 
connects to its nearest BS under the proposed UE sampling 
method with thinning probability p is 

-l 



Pc(a,/3,p) 



1 



which for a = 
P C (4,/3,P) 



p^ 
4 simplifies to 

p^ (tt/2 



1 



-s 1 



-du 



u 2 
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This result directly follows from Theorem 2 of Q with a 
slight modification to incorporate conditional thinning of the 
point process <£>. The coverage probability is first conditioned 
on the distance of a typical UE to its serving BS, which can 
be computed by the null probability of a PPR Conditioned on 
this distance, say u, the interference field defined over R 2 n 
B(Q, u) c is a PPP with density pX, where B(0,u) is a ball 
with radius u centered at and p appears due to conditional 
thinning. This is the step where the proof of the current Lemma 
differs from that of Theorem 2 of [ 1 1 . The remaining proof, 
which mainly involves the derivation of the Laplace transform 
of interference, remains the same. 

Remark 2 (Scale invariance and effect of p on P c ). For any 

given value of thinning probability p, the coverage probability 
is independent of the density of the BSs in an interference 
limited cellular network. This scale invariance result is a 
generalization of a similar observation reported in ^ for the 
uniform distribution of the UEs, which is a special case of the 
proposed model and corresponds to p = 1. As expected, the 
coverage probability is a monotonically decreasing function 
of p. This is consistent with the observations reported in / [70y , 
where a similar parametrization was used to model load on 
various classes of BSs in a heterogeneous cellular network. 




Fig. 3. A realization of the proposed non-uniform UE distribution model. 



IV. Non-Uniform UE Distribution 

In this section, we show that the framework developed in the 
previous section can be used as an accurate analytical genera- 
tive model for the non-uniform UE distributions where the UEs 
are clustered around their serving BSs. We first explain the 
simulation setup for this non-uniform UE distribution model 



and then describe a subtle difference in this model and the UE 
sampling framework developed in the previous section. 

To simulate the non-uniform UE distribution model, start 
with a realization of a PPP $ with density A and form the 
Voronoi tessellation. Recall that in case of a cellular network, 
the Voronoi cell of each point denotes its coverage region 
in the nearest neighbor connectivity model. Distribute N u 
users uniformly in each Voronoi cell. For concreteness, we 
assume that N u is the same for all the BSs and equal to the 
number of resource blocks, which models a full-buffer system. 
This assumption can be relaxed under certain conditions as 
discussed in flO) . Until this point, the UE distribution is fairly 
uniform although there is a subtle difference in this model 
and the way typical user is usually defined to be uniformly 
distributed over W 2 . We will remark on this difference later in 
this section. To induce dependence in the BS and UE point 
processes, we again use the the thinning idea and retain points 
of the realization of $ independently with probability p and 
remove the rest. The UEs corresponding to the points that 
are removed are also removed. The thinned version of the 
point process <£> with density \p models the BS locations. For 
the new coverage areas defined by the thinned point process, 
the remaining UEs are biased towards the cell interior as 
shown in Fig. [3] A favorable characteristic of this model 
is the probabilistic attraction introduced in the UE and BS 
point processes without inducing any geometric constraints. 
The coverage probability can now be numerically evaluated 
by averaging over these UE locations. Note that the result of 
Lemma [2] is not exact for this simulation model, as remarked 
below. 

Remark 3 (UE uniformly distributed in IR 2 vs in randomly 
chosen Voronoi Cell). There is a subtle difference in perform- 
ing downlink analysis at a typical UE uniformly distributed in 
M 2 and uniformly distributed in a randomly chosen Voronoi 
cell. The former corresponds to the analytical model discussed 
in the previous section for p = 1 and the latter corresponds to 
the simulation model discussed in this section. The difference 
is induced by the structure of Poisson Voronoi tessellation and 
can be understood by a simple fact that a point uniformly 
distributed in IR 2 is more likely to fall in a bigger Voronoi 
cell, whereas there is no such bias when we choose a point 
uniformly distributed in a randomly chosen Voronoi cell. 

We now compare the coverage probability of this non- 
uniform UE distribution model with the analytical expression 
derived in Lemma [2] in Fig. [4] We first note that the plots are 
surprisingly close and the difference highlighted in Remark [3] 
has a negligible impact on the coverage probability. Thus, the 
analytical model based on conditional thinning can be used 
as an accurate generative model to study coverage probability 
for this non-uniform UE distribution model described in this 
section. Further, we compare the coverage probability with the 
baseline model, where the UEs are distributed uniformly over 
K 2 independent of the BS point process. We note that the 
difference in coverage predictions is significant even for high 
values of p. This clearly highlights the importance of accurate 
UE distribution models in the performance analysis of cellular 
networks. 




Thinning Probability (p) 

Fig. 4. Comparison of the coverage probability of the proposed non-uniform 
UE distribution model with the analytical expression derived in Lemma[2]and 
the baseline model assuming uniform UE distribution (a = 4). 



V. Conclusion 

In this letter, we addressed the problem of incorporating 
non-uniform UE distributions in the random spatial models 
for cellular networks. This work has numerous extensions, 
both in system modeling and analysis. In system modeling, 
it is important to develop empirical UE distribution models, 
especially for heterogeneous networks, where the UEs are 
more likely to be around low power nodes. In analysis, it 
is important to develop tractable modeling tools capable of 
handling correlations in the UE and BS point processes. 
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